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Abstract 

We investigate the persistence probability in the Voter model for 
dimensions d>2. This is achieved by mapping the Voter model onto 
a continuum reaction-diffusion system. Using path integral meth- 
ods, we compute the persistence probability r(q,t), where q is the 
number of "opinions" in the original Voter model. We find r(q,t) ~ 
exp[-/ 2 (g)(mt) 2 ] ind = 2; r(q,t) ~ exp[- f d {q)t {d ~ 2)l2 \ for 2 < d < 4; 
r(q,t) ~ exp[-/ 4 (g)t/lni] in d = 4; and r(q,t) ~ exp[-f d (q)t] for 
d > 4. The results of our analysis are checked by Monte Carlo simu- 
lations. 

The Voter model is a simple stochastic model which exhibits interesting di- 
mension dependent properties 0. Whilst, in one dimension, it is equivalent 
to the Glauber-Ising model at zero temperature, its properties depart from 
this model in higher dimensions. On each site of a d— dimensional lattice, 
opinions of a voter, or values of a spin a = 1,2, ... ,q, are initially distributed 
randomly. Between t and t + dt a site is picked at random. The voter on this 
site takes the opinion of one of its 2d neighbours, also chosen at random. This 
model is equivalent to a system of coalescing random walkers II] . Therefore, 
due to the recurrence properties of random walks, the interactions between 



walkers are strong in low dimensions (d < 2), and progressively decrease in 
higher dimensions. As a consequence, for d < 2 a, consensus of opinion is 
approached as t — > oo, or in other words, the system coarsens. This con- 
trasts with the behaviour in higher dimensions where disagreements persist 
indefinitely 

In this letter we compute the persistence probability r(q, t), i.e. the prob- 
ability that a given voter never changed his opinion up to time t, in arbitrary 
dimension d > 2, using field theoretical methods. Persistence probabilities 
have been investigated for a number of coarsening systems, including the 
Glauber-Ising chain [Q, [3], [|], exhibiting in all cases an algebraic decay with 
time. By contrast, it was found numerically that the persistence probability 
of the 2c? Voter model with q = 2 had the behaviour exp[— const. (In t) 2 } [[J. 
This change in behaviour with dimension motivated the present work. Field 
theoretical methods are well adapted for the study of the role of dimension 
in simple physical systems, hence their use is natural in the present context. 

The first step for the study of persistence consists in mapping the Voter 
model onto its dual particle system, i.e. a system of particles on a d— dimensional 
lattice, with the reaction A + A —>■ A, fermionic occupation rules and parti- 
cle injection at the origin. As soon as one particle leaves that site another 
particle is immediately injected there in its place. This reaction-diffusion 
model is hereafter referred to as model I. Defining Q(m,t) as the prob- 
ability that, after a time t there are m particles in the system, one has 
r(q,t) = Era=i Q(jn,t) q l ~ m 0, f|. After a time t the injected particles will 
have penetrated a distance of order t 1 / 2 from the source. At greater distances 
there are almost no particles present, whereas behind the front a steady state 
is reached. Hence in order to calculate the probability r(q, t), one can proceed 
as follows: (i) Compute 

N 

s(q,L)=Y,P\rn,L)q 1 - m , (1) 

m=l 

where P I (m, L) is the steady state probability of finding m particles within 
a distance L of the point of injection (for t 1 / 2 ^> L), and where N = 0(L d ) 
is the number of sites within that radius, (ii) Now compute the persistence 
probability using finite size scaling: r(q,t) ~ s(q, L ~ t 1 / 2 ). 

Unfortunately a computation of the P 1 (m, L) is still too difficult. There- 
fore we replace model I with a slightly modified variant, which we will call 
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model II. We now have multiple (bosonic) site occupancy with a finite re- 
action rate A, and particle input at the origin with a rate J. However the 
relevant properties of both models will turn out to be dominated by their 
(low density) behaviour far from the injection point. This means that, in 
the large L limit, the distributions P 1 and P 11 actually become very similar, 
hence model II can be used for the persistence calculation. Evidence for this 
claim will be given later. 

The mean field equation for model II, with particle density p(x, t), is given 
by (for details see f| [7j) 

d tP = DV 2 p- Xp 2 + J6 d (x), (2) 

where D is the diffusion constant. Well behind the diffusive front, but far 
from the injection point, we can insert a steady-state power law ansatz into 
(D, giving 

(2(A-d)D/(\x 2 ), (d<4) 
p{x) ~ | 2D/(Xx 2 \nx), (d = 4) (3) 
[ const. x 2 ~ d , (d > 4). 

Note that for d < 4, the asymptotic form of p(x) does not depend on J, an 
indication of the similarity between models I and II. The change in behaviour 
at d = 4 is due to the role of the reaction term in (0). For large x, this reaction 
term is marginal in d = 4, whereas for d > 4, it can be asymptotically ignored, 
the density instead decaying in accordance with Laplace's equation. Hence 
d = 4 is, in a sense, a critical dimension. Nevertheless the behaviour of the 
average density around d = 4 can be described adequately within mean field 
theory. However in low spatial dimensions (d < d c — 2), the mean field 
equation (Q) is no longer correct, since it does not include fluctuation effects. 
We therefore map model II to a field theory which systematically includes 
these fluctuations, and permits an analysis of the probability distribution 
P u (m, L) in arbitrary dimension. Following the work of || M, and using the 
formalism of HlQ[ , the master equation for model II may be rewritten as a 



second quantized Hamiltonian 

H = -Dj2{al(aj - a,)} - A£{(1 - aMd 2 } - J{d\ - 1}. (4) 

n.n. i 

The first term describes diffusion of the A particles, the second term de- 
scribes their coagulation reaction, and the final term results from stochastic 
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(Poissonian) particle injection at the origin (site 1). Defining P{n\, . . . ,un) 
to be the probability that the system is in a configuration {m, . . . , un} {ni is 
the occupation number for site i), then the steady state probability P u (m, L) 
of finding m particles in the system is 



P u (m, L)=Y1 P ( ni ' •••■> u n) 5(52 n i> m )- 



(5) 



{n} 



In the second quantized formalism, where is the system state ket, this 
becomes 



P u (m,L) = (0|exp($» <$($>,-, m) |tf> (6) 

(0\a^...a n N N \V). 



- E 

to' ^ 



3 

to! 



ni! . . . n^! 



n,A — m 



Identifying this last expression as a multinomial expansion yields 

P ll (m,L) = (0\^-(52a l ) m m- (7) 
Hence the generating function for the probabilities P n (m, L) is given by 

oo 

G(y) = £ y m P n (m,L) = (0| «p(y5»l*> = (i «p[(l/ - 1) ( 8 ) 

m=0 i i 

with the definition (| = (0| exp(J2i Performing the mapping to a path 
integral || |9|, yields the continuum action 



S = J d d x J dt (^d(d t - V 2 )a + Aaa 2 + AaV - Ja5 d (x] 



(9) 



where we have absorbed the diffusion constant D into a rescaling of time and 
of the couplings A and J. Furthermore the continuum limit of (§) is 

G(y) = (exp [(</ -I) J d d x a(x)] ^ , (10) 

with the average performed with respect to the action S. From (|8]) we have 

1 d ri 



P n {m, L) 



ml dy 7 



My) 



(11) 



y=0 



which leads to the continuum expression 
1 



P u (m,L) 



ml 



d d Xi a(xi) 



- / d x a(x) 



and hence, from ([I]), 



s(q,L) ~ G(q r ) = (exp 



d x a(x] 



;i2) 



(13) 



where we have used the equivalence between models I and II in the large L 
limit, mentioned previously. 

The reaction rate A has naive scaling dimension k 2 ~ d . This predicts a 
critical dimension d c — 2 at or below which we must renormalize the the- 
ory. A very similar case has been studied in [IT]. Hence the dimensionless 
renormalized coupling evaluated at the normalization point s = k 2 is given 
by g R = n~ e \R, where e = 2 — d and 

X R = \(l + \I d (K))-\ /„(«) = T' 2 



d d p 



(2n) d s+p 2 



(14) 



The beta function is then given by f3(gn) = nd K gR = g^/27T in d = 2. This 
is exact to all orders in g^. The field (a) now satisfies the renormalization 
group equation 



(a)(x,g R ,K) 



0. 



(15) 



The solution of this equation is given by 

(a)(x,g R , «) = (Kx)' d (a)(K'' 1 ,g R , k) (d<2), (16) 

where the running coupling ^ acquires the universal asymptotic form 

2tt 



lnx 



l + OfOna;)- 1 ] (d = 2). 



(17) 



The leading order (tree) diagrams for the density are exactly equivalent to 
the mean field rate equations |10| . Hence inserting the mean field expression 
(H) into the RG equation (|16|) yields the exact leading asymptotic expression 



(a(x)) 



2 lnx 



7rar 



(l + 0[(\nx)- 1 }) (d = 2). 
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(18) 



We note that this renormalization procedure is only needed at or below the 
upper critical dimension. For d > 2 the mean field results in (||) become 
qualitatively correct. 

The next step is to turn the expression ( |T3"D for the persistence probability 
into a cumulant expansion (see |12| for a similar case), 



s(q,L) 



exp 



1 -q- 1 ) J d d x a(x) 



(19) 



exp 



-(1-g- 1 ) J d d x (a(x)) + ^(1 - q- 1 ) 2 J d d x J d d x' (a(x)a(x')) c 



We now consider the behaviour of s(q, L) in different dimensions. 

• d = 2. We use fll8[) to derive the leading behaviour of the first term in the 

cumulant expansion. Hence for a system of radius L, 



d 2 x (a(x)) 



d x 



2 lnx 



7TX Z 



2(\nLf 



(20) 



where the next order term is suppressed by (InL) . This result implies that 
almost all of the particles in the system are situated very far from the point of 
injection. This confirms that the behaviour of the model close to the injection 
point is unimportant when considering the large L limit. Substituting ( p0|) 
into (|I9| ) gives 



s(q,L) - exp [-h 2 (q)(\nL) 2 + h {2) (q) \nL + 
h 2 (q) = 2(1 - q- 1 ) + 4C 2 (1 - q' 1 ) 2 + ■■■ , 



(21) 
(22) 



where we assume that all the higher order cumulants also have leading order 
contributions of (9(lnL) 2 . Simulations of the reaction-diffusion system have 
confirmed this expectation (see below). Hence for the Voter model, using 
finite size scaling L ~ t 1 ! 2 in (|2T|), we have 



r(q,t) ~ exp [-/ 2 (g)(lnt) 2 + jf\q) lnt + 
/ 2 (g) = (l/2)(l-g- 1 )+C 2 (l-g- 1 ) 2 + - 



(d = 2) (23) 
(24) 



This analysis shows that the unusual logarithmic form of this decay is due 
to the fact that d = 2 is a critical dimension. 
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• 2 < d < 4. Using the mean field expressions ([3]) yields the leading behaviour 
of the first term in the cumulant expansion 

Ad i , w f L d - 2 /\, (2<d<4) 
fx (a(x)) oc | LV fa L)i \ d = ^ (25) 

Hence almost all of the particles are again situated far from the injection 
site. Furthermore, dimensional analysis can be combined with an analysis of 
the diagrams for the higher order cumulants to show that each term in the 
cumulant expansion scales identically to (|25|) . We assume that a logarithm 
is present in each term of the cumulant expansion for d = 4, although this 
analysis only explicitly shows the presence of a logarithm in the first term. 
Using finite size scaling, we then have 

-f d (q)t (d - 2)/2 + ■ ■ ] , (2<rf<4) m) 
-f A (q)t/\nt + ■■■}, (d = 4) 

g- 1 ) + Q(l-g- 1 ) 2 + ---. (27) 

Note that the cumulant amplitudes B d , Cj, ■ ■ ■ are no longer universal. 

• d > 4. We again use (§) for the first term in the cumulant expansion, 

d d x (a(x)) oc L 2 . (28) 

For d > 4 the reaction process is asymptotically irrelevant since the particles 
no longer "see" each other. Since the (non-Poissonian) reaction process is not 
then present, the diagrams for the higher order cumulants actually vanish. 
Hence we have 

r(g,t)~exp[-/ rf (g)t + ...] (d > 4) (29) 
f d (q) = B d (l-q- 1 ). (30) 

In summary, we find that, at large times and in any dimension, r(q,t) ~ 
exp[— f(q)(m(t))], where (m(t)) is the average number of particles in the 
dual reaction-diffusion models at time t (see also ||13|| ). Note that in d = 1, 
(m(t)) grows logarithmically in time, in accord with the algebraic fall off of 
the persistence probability. 

Let us finally report on simulations. We measured the persistence prob- 
ability r(q, t) for the Voter model, as a function of time, with q = 2, . . . , oo. 
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The results for d = 2 on a system of size (1000) 2 are shown in figure 1, 
where a decay r(q,t) ~ exp[— f2(q)(\nt) 2 ] is clearly seen, in agreement with 
both our predictions and the simulations for q — 2 in JJ. A plot of f^il) 
extracted numerically is given in figure 2. It shows qualitative agreement 
with the prediction given in ([24]) . In particular for increasing values of q, 
fz(q) progressively departs from the first predicted term (dashed line on the 
figure). 

Simulations performed on the 3c? Voter model of size (100) 3 do not show 
the expected decay r(q,t) ~ exp[— f2,(q)t 1 / 2 ] at large times. Instead, the 
numerical results indicate an effective exponent for t slightly larger than 1/2. 
This is due to large subleading corrections which make extraction of the 
leading behaviour difficult. 

We now turn to simulations of the dual model I. In d = 2 we measured 
(to) = J d d x{a(x)), the average number of particles in the stationary state, 
for system sizes L = 10, ... , 80, averaged over 2 x 10 6 MC steps (figure 3). 
The data are consistent with (m) ~ 2(lnL) 2 plus corrections of (9(lnL), in 
agreement with (p0|). This plot also provides evidence for the equivalence 
between models I and II, since our simulations are for model I whereas the 
theory was developed for model II. 

We also examined the behaviour in d = 2 of the first three higher order 
cumulants for the a(x) field (integrated over the system size as in (^)). They 
all scale with L in the same way as the first cumulant / d d x{a(x)) = (to). 
This confirms that, in principle, each term in the cumulant expansion is 
equally "relevant" for the persistence calculation. Nevertheless the numerical 
prefactors for each term appear to decrease fairly rapidly indicating that the 
first few terms are the most significant. 

Finally, in d = 3 we measured (to) as a function of L, for L — 10, ... , 40. 
For large sizes, instead of the prediction (m) ~ L (see (p5|)), the data show 
an apparent exponent of about 1.2, indicating again the presence of large 
subleading corrections. 
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Figure 1: lnr(g,t) plotted against (hit) 2 for q = 2, 3, 4, 5, 6, 7, 8, 10, 14, 20, 
40, 10 6 (from top to bottom), in the 2d Voter model (system size (1000) 2 ). 
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Figure 2: Plot of the numerical data for f2{q), extracted from figure 1, versus 
(1— q^ 1 ) (dots). Also illustrated is the lowest order theoretical approximation 
t° J2(q) (dashed line). 
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Figure 3: Plot of (m) versus InL for system sizes L — 10, ... , 80, averaged 
over 2 x 10 6 MC steps (for model I). Solid line is a quadratic best fit to the 
data, where the coefficient of the quadratic term is found to be close to 1.9. 
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